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NONLINEAR DYNAMICS OF GLOB AT, ATMOSPHERIC AND 
EARTH SYSTEM PROCESSES 


During the grant period we continued our ongoing studies aimed at enhancing our 
understanding of the operation of the atmosphere as a complex nonlinear system interacting with 
the hydrosphere, biosphere, and cryosphere in response to external radiative forcing. Five papers 
were completed with support from the grant (NAG 5-2316), representing contributions in three 
main areas of study: I) theoretical studies of the interactive atmospheric response to changed 
biospheric boundary conditions measurable from satellites, II) statistical-observational studies of 
global-scale temperature variability on interannual to century time scales, and III) dynamics of 
long-term earth system changes associated with ice sheet surges. These works involve the efforts 
of T. Zhang, M. E. Mann, M. Verbitsky, and the principal investigator. A listing of the papers is 
followed by the title pages and abstracts of the listed papers. 
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Sensitivity Properties of a Biosphere Model Based on BATS 
and a Statistical-Dynamical Climate Model 

Tajping Zhang 

Department of Geology and Geophysics, Yale University, New Haven . Connecticut 
(Manuscript received 23 February 1993, in final form 29 October 1993) 

ABSTRACT 

A biosphere model based on the Biosphere-Atmosphere Transfer Scheme (BATS) and the SalUman-Vemekar 
(SV) statistical-dynamical climate model is developed. Some equations of BATS are adopted either intact or 
with modifications, some are conceptually modified, and still others are replaced with equations of the SV 
model. 

The model is designed so that it can be run independently as long as the parameters related to the physiology 
and physiognomy of the vegetation, the atmospheric conditions, solar radiation, and soil conditions are given. 
With this stand-alone biosphere model, a series of sensitivity investigations, particularly the model sensitivity 
to fractional area of vegetation cover, soil surface water availability, and solar radiation for different types of 
vegetation, were conducted as a first step. These numerical experiments indicate that the presence of a vegetation 
cover greatly enhances the exchanges of momentum, water vapor, and energy between the atmosphere and the 
surface of the earth. An interesting result is that a dense and thick vegetation cover tends to serve as an environment 
conditioner or, more specifically, a thermostat and a humidistat, since the soil surface temperature, foliage 
temperature, and temperature and vapor pressure of air within the foliage are practically insensitive to variation 
of soil surface water availability and even solar radiation within a wide range. An attempt is also made to 
simulate the gradual deterioration of environment accompanying gradual degradation of a tropical forest to 
grasslands. Comparison with field data shows that this model can realistically simulate the land surface processes 
involving biospheric variations. 


1. Introduction 

For a realistic simulation of the climate, the atmo- 
sphere, hydrosphere, cryosphere, lithosphere, bio- 
sphere, and their interactions with one another should 
be described simultaneously on the basis of physics. 
Great achievements have been made in the develop- 
ment and application of general circulation models 
(GCMs), but many inadequacies remain because of 
an incomplete understanding of the physics involved 
and inappropriate parameterizations of some subgrid 
processes. The land surface processes, or the interac- 
tions between the atmosphere and biosphere, are just 
such processes to which the climate modeling com- 
munity had not paid much attention until several years 
ago. There have now been some studies and observa- 
tions that indicate the existence and the spatial and 
temporal variation of vegetation on the continent have 
consequential effects on the maintenance and change 
of the global climate (Shukla and Mintz 1982; Mintz 
1984; Henderson-Sellers 1990; Shukla et al. 1990; 
Nobreetal. 1991; Gash 1991). 
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The vegetation on the land surface, by altering the 
surface albedo, the drag coefficient, and the partitioning 
of available energy into sensible and latent heat fluxes, 
greatly complicates the surface processes, particularly 
the exchange of energy, mass, and momentum between 
the surface and the atmosphere. The maintenance of 
soil moisture is also closely related to the existence of 
vegetation cover. Attempts have been made to describe 
and parameterize the surface processes, but even the 
most sophisticated models give only a first approximate 
estimate. In addition, the requirement for satisfactory 
modeling of the surface processes at much higher res- 
olutions than is feasible for GCMs gives rise to added 
difficulties in incorporating the biosphere model into 
the GCMs. 

The state-of-the-art representatives of biosphere 
models for GCMs are the Simple Biosphere Model 
(SiB) of Sellers and Mintz ( 1986), the Biosphere-At- 
mosphere Transfer Scheme (BATS) of Dickinson et 
al. ( 1986), and the Bare Essentials Model of Pitman 
( 1988), which have been incorporated into GCMs of 
the Goddard Laboratory for Atmospheric Sciences, the 
National Center for Atmosphere Research, and the 
Canadian Climate Centre, respectively. These sub- 
models, though elaborate compared with previous 
ones, are conceptually much simpler than reality and 
are still so-called “big-leaf” models, a concept first in- 
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Heinrich-type glacial surges in a low-order dynamical climate 
model 

Mikhail Verbitsky, Barry Saltzman 
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Received 12 July 1993/ Accepted 15 October 1993 


Abstract. Recent studies suggest the occurrence of spo- 
radic episodes during which the ice streams that dis- 
charge ice sheets become enormously active, produc- 
ing large numbers of icebergs (reflected in North At- 
lantic sea cores as “Heinrich events") and possibly 
causing the partial collapse of the ice sheets. To simu- 
late the mechanism of internal thermo-hvdrodynamical 
instability implied by such behavior in the context of a 
more general paleoclimate dynamics model (PDM), we 
introduce a new sliding-catastrophe function that can 
account for ice-sheet surges. In particular, using simple 
scaling estimates derived from the equations of motion 
and thermo-conductivity for ice flow, we express this 
function in terms of the thickness, density, viscosity, 
heat-capacity, and heat-conductivity of ice. Analysis of 
the properties of this function suggests that these Hein- 
rich-type instability events might be of three possible 
kinds; the first type of event occurs in periods of glacial 
maximum when temperature conditions on the ice sur- 
face are extremely cold, but internal friction within 
bottom boundary layer is also at its maximum and is 
strong enough to melt ice and cause its surge. The sec- 
ond type of event may happen during an interglacial, 
when the ice thickness is small but relatively warm cli- 
matic conditions on the upper surface of ice can be eas- 
ily advected with the flow- of ice to the bottom where 
even a small additional heating due to friction may 
cause melting. The third and. perhaps, most interesting 
type of event is one that may occur during ice sheet 
growth; in this period particles of ice reaching the bot- 
tom still “remember" the warm temperature condi- 
tions of the previous interglacial and additional heating 
due to increasing friction associated with the growing 
ice sheet may again cause melting. To the extent that 
the upper glacier surface temperature depends on at- 
mospheric carbon dioxide concentration, this third 
case introduces the interesting possibility that earlier 
C0 2 concentrations may be as important for the pres- 
ent-day climate as its current value. We present results 
of numerical experiments demonstrating how these 
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three kinds of instability can originate and interact 
with other components of the global climate system to 
produce variations of the Heinrich-event type. In parti- 
cular, according to our model the climate system seems 
more vulnerable to surges during the penultimate in- 
terglacial period than in the present one, which may 
contribute to an explanation of the recent results of the 
Greenland Ice Core Project. 


Introduction 

Recent geological studies (e.g., Andrews and Tedesco 

1992) relate pulses of detrital carbonate-rich sediments 
in cores from the Labrador Sea to episodes in the late 
Pleistocene history of the Laurentide ice sheet when 
ice streams from the Hudson Strait become active 
enough to be able to deliver sediments onto the deep- 
sea plain. These incidents correlate highly with “Hein- 
rich events" in North Atlantic sedimentary cores 
(Heinrich 1988; Broecker et al. 1992; Grousset et al. 

1993) . Another, theoretical, study (MacAyeal 1992) 
suggests that the West Antarctic ice sheet experienced 
sporadic collapses throughout its last-million-year's 
evolution. Although the first group of events has a pe- 
riodicity on the order of 10 ky and the second group 
(according to the numerical experiments) took place 
only a few times during the last million years, we shall 
treat them both as Heinrich-type events because it is 
not unlikely that both of these types of episodic behav- 
ior have the same physical explanation, i.e., internal 
thermo-hvdrodynamical instability inherent to the 
large ice sheets due to basal processes (e.g., Paterson 
1981; Hughes 1992). 

The aim of our study here is (1) to describe this 
kind of instability on the basis of simple scaling esti- 
mates derived from the motion and thermo-conductivi- 
tv equations of large ice masses. (2) to enlarge the phy- 
sical content of our previous paleoclimate dynamics 
model (PDM) described in Saltzman and Verbitsky, 
i.e., S-Verb 1993. by including the possibility for this 
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Global-scale modes of surface temperature variability 
on interannual to century timescales 


Michael E. Mann and Jeffrey Park 


Department of Geology and Geophysics. Yale University, New Haven, Connecticut 

Abstract. Using 100 years of global temperature anomaly data, we have performed 
a singular value decomposition of temperature variations in narrow frequency 
bands to isolate coherent spatio-temporal u modes ,? of global climate variability. 
Statistical significance is determined from confidence limits obtained by Monte 
Carlo simulations. Secular variance is dominated by a globally coherent trend, 
with nearly all grid points warming in phase at varying amplitude. A smaller, 
but significant, share of the secular variance corresponds to a pattern dominated 
by warming and subsequent cooling in the high latitude North Atlantic with a 
roughly centennial timescale. Spatial patterns associated with significant peaks in 
variance within a broad period range from 2.8 to 5.7 years exhibit characteristic 
El Nino-Southern Oscillation (ENSO) patterns. A recent transition to a regime 
of higher ENSO frequency is suggested by our analysis. An interdecadal mode in 
the 15-to-18 years period range appears to represent long-term ENSO variability. 
This mode has a sizeable projection onto global-average temperature, and accounts 
for much of the anomalous global warmth of the 1980s. A quasi-biennial mode 
centered near 2.2-years period and a mode centered at 7-to-8 years period both 
exhibit predominantly a North Atlantic Oscillation (NAO) temperature pattern. 

A potentially significant MecadaU mode centered on ll-to-12 years period also 
exhibits an NAO temperature pattern and may be modulated by the century-scale 
North Atlantic variability. 



1. Introduction 

In the face of possible anthropogenic effects on global 
climate, there is a need to characterize better the na- 
ture of historical climate variability. Aside from possi- 
ble episodic volcanic forcing [e.g., Bradley and Jones , 
1992], organized interannual climate variability seems 
to be associated largely with the El Nino-Southern 
Oscillation (ENSO) [e.g., Philander , 1990], the quasi- 
biennial oscillation (QBO) [e.g., Trenberih and Shin , 
1984], and regional climate processes, such as the North 
Atlantic Oscillation (NAO) [e.g., Bjerknes , 1964; van 
Loon and Rogers , 1978; Rogers , 1984; Lamb and Pep- 
pier , 1987] which may have hemisphere- or even global- 
scale influence. Decadal and interdecadal variability is 
less well understood and has variously been attributed 
to internal variability in the thermohaline circulation 
of the oceans [ Levitus , 1989; Ghil and Vautard , 1991; 
Stocker and Mysak , 1992; Read and Gould , 1992], cou- 
pled cryosphere-ocean-atmosphere cycles [Mysak and 
Power , 1992], and external astronomical forcing due 
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to the 18.6-year soli-lunar tide [ Cume and O'Brien , 
1992; Mitra et ai, 1991], the ~ 11-year solar cycle [Lab- 
itzke and van Loon . 1988: Mitra et a/., 1991; Cume and 
O'Brien, 1992] and its 22-year subharmonic or '‘Hale” 
cycle [Vines, 1986]. On century and longer timescales, 
thermohaline and cryosphere variability [Stocker and 
Mysak , 1992], solar forcing [Friis- Christensen and Las- 
sen , 1991], and anthropogenic factors [e.g., Interna- 
tional Panel on Climate Control (IPCC), 1990] have 
been argued to be important. 

While there is scant evidence for truly periodic histor- 
ical interannual climate signals (that is, phase-coherent 
sinusoidal variations) there is nonetheless strong evi- 
dence that certain climatic processes are quasi-periodic. 
Such processes exhibit well-defined timescales of oscil- 
lation but are frequency or amplitude modulated on 
a longer timescale, often in an unpredictable manner. 
Thus spectral energy is spread out in a finite band about 
some dominant frequency. For example, ENSO variabil- 
ity appears to concentrate into distinct low-frequency 
(4-6 year) and high-frequency (2-3 year) bands [e.g., 
Barnett , 1991; Keppenne and Ghil, 1992; Dickey et al . , 
1992; Ropelewski et ai, 1992]. Quasi-periodic modes of 
climate variability, then, should be identifiable as sta- 
tistically significant, relatively narrowband peaks in the 
spectra of global climate records. Any global modes of 
climate variability should also be associated with spa- 
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GLOBALLY CORRELATED VARIABILITY IN SURFACE TEMPERATURES 

Michael E. Mann and Jeffrey Park 

Dept, of Geology and Geophysics j - 

Yale University 
New Haven, Connecticut 


1. INTRODUCTION 

Several workers (eg., Folland et al 1984; Ghil and Vau- 
tard, 1991) have examined records of globally- and 
hemispherically-averaged land air and/or sea surface 
temperature anomalies for signals on interannual to 
century time scales. Mann and Park (1993) (herafter 
MP93) demonstrated that interdecadal temperature 
variability exhibits spatial variability that leads to sig- 
nificant cancelation in a global average. An analysis 
based only on modes which have large projections onto 
global-average temperature thus gives an incomplete, 
and perhaps misleading, picture of global-scale vari- 
ability. We seek to isolate the spatial modes of global 
temperature variability on interannual to century-long 
time scales. This approach allows us to examine cli- 
mate processes associated with changes in atmospheric 
circulation that lead primarily to a redistribution of 
heat over the surface of the earth, as well as those 
related to changes in the overall surface heat budget. 

Quasi- periodic variability, such as the El-Nino 
/Southern Oscillation (ENSO), can be distinguished 
by correlated patterns of variability, or teleconncc- 
tions, between widely spaced portions of the Earth's 
surface. Our data analysis seeks such patterns that arc 
statistically distinguished from a locally white-noise 
background. 

2. METHOD 

We use gridded land air and sea surface temperature 
anomalies distributed on a 5° longitude by 5° latitude 
grid (see Jones and BrifFa, 1992) spanning N 1~00 
months from 1891 through 1990, keeping only a subset 
(M = 449 grid points) with nearly continuous monthly 
sampling. Potential biases in the data set are discussed 
elsewhere (eg., Jones 1986). Because the data samples 
the Earth's surface unevenly, poorly sampled regions 
may be weakly represented in the analysis, while well- 
sampled regions may exert influence out of proportion 
with its true place in the global climate system. Few 
historical climate datasets of significant length, how- 
ever, are free from this shortcoming. 

We demean and normalize the rnth grid point time 
series x^ (m) (n = 1 , . . . , A months) by its monthly 
standard deviation to assure that regions with 
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high variance (e.g., continent centers) do not in- 
ordinately influence our analysis. To study quasi- 
periodic phenomena, we transform from the time 
to the spectral domain using multitaper spectral 
analysis (Thompson, 1982). Given the set of nor- 

1 ■ (rn) _ ' ( m ) /^( m ) we cai- 

malized time series, Xt» — x n i a 

culate the tapered Fourier transforms Y k m (/) = 

where At = 1 month is the 

sampling interval and {u)^}n=i ‘ s member 

in an orthogonal sequence of Slepian tapers. A set 
of k tapers that have power concentrated over a half- 
bandwidth of p/r are referred to as p7r-prolate Slepian 
tapers, where fn = (N At)~ l is the Rayleigh fre- 
quency. The choice of K represents a compromise be- 
tween the variance and resolution of the spectral repre- 
sentation, as the tapers average over a half-bandwidth 
of pf R and only the first 2 p - 1 tapers are usefully 
resistant to spectral leakage. 

Assuming a locally white spectral background, the 
multitaper technique produces A tapered Fourier 
transforms at a given frequency / for each of the M 
grid points, from which we construct the matrix 


vV *> 

y a (,) • 

. Vk (1) ‘ 

Yi™ 

vY 2) • 

. w (2) 

y/(M) 

vY M > . 

.. y k ' m) . 


We decompose A with a complex-valued SVD into A 
orthonormal left-eigenvectors u fc , representing empiri- 
cal orthogonal functions (EOFs) in the spatial domain, 
and K orthonormal right-eigenvectors v fc , representing 
EOFs in the spectral domain. The v* are complex- 
valued A'-vectors, and can be inverted to obtain the 
slowly- vary ing envelope of the kth mode. The 
complex-valued M-vectors which determine the am- 
plitude and relative phase of temperature variations 
among grid points, The singular value X k scales the 
amplitude of the k th mode. We order the singular 

values Ai(/)>A 2 (/)^---^^h(/) — 

Park (1992) shows how the slowly-varying envelope 
A(<) of a quasi-periodic signal x(t) = ^{A(t)e } 

can be estimated from a set of eigenspectra Y k (f), k = 
1 AN In the multivariate case, the k th time-domain 
signal is reconstructed from v fc (/), with some addi- 
tional constraints. We minimize the first-derivative 
of the discrete representation of the envelope A(t) 
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Decadal-to- Centennial- Scale Climate Variability: 
Insights into The Rise and Fall of the Great Salt Lake 

Michael E. Mann 

Department of Geology and Geophysics, Yale University, New Haven, Connecticut 

Upmanu Lall 

Utah Water Research Laboratory, Utah State University, Logan, Utah 

Barry Saltzman 

Department of Geology and Geophysics, Yale University, New Haven, Connecticut 

Abstract. We demonstrate connections between 
decadal and secular time scale global climatic varia- 
tions, and historical variations in the volume of the 
Great Salt Lake. The decadal variations correspond 
to a low-frequency north-south shifting of storm tracks 
which influence winter precipitation and explain nearly 
18% of the interannual and longer-term variance in the 
record of monthly volume change. The secular trend 
accounts for a more modest ~1.5% of the variance. 


Introduction 

The prospect of anthropogenic climate change has 
renewed interest in the study of low frequency vari- 
ability in hydroclimatic processes [Aguado et al 1992; 
Rasmusson and Arkin, 1993; Lins and Michaels , 1994]. 
Such concerns are especially pertinent with respect to 
water availability in arid regions. Permanent closed 
basin lakes are unique in that they reflect a precar- 
ious long term balance of precipitation and evapora- 
tion in an arid region, and amplify persistent imbal- 
ances. Several studies have established linkages between 
regional hydrological patterns and large-scale atmo- 
spheric circulation [Rasmusson and Arkin , 1993; Lins , 
1985; Cayan and Peterson , 1989; Ely et al 1993; 1994]. 
Here, we demonstrate the connection between large- 
scale climatic signals and the fluctuations in the vol- 
ume of the closed-basin Great Salt Lake (GSL). The 
GSL (Figure la) provides a particularly useful hydro- 
logical case study because of its dramatic historical 
fluctuations (Figure lb), emphasizing lower frequencies 
through a spatial and temporal integration of climatic 
processes. The monthly volume change (Al/-Figure lc) 
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